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Abstract 

Recently it has been shown that a large fraction of the dwarf satellite galaxies orbiting the An¬ 
dromeda galaxy are surprisingly aligned in a thin, extended and kinematically coherent planar struc¬ 
ture. The presence of such a structure seems to challenge the current Cold Dark Matter paradigm of 
structure formation, which predicts a more uniform distribution of satellites around central objects. 

We show that it is possible to obtain a thin, extended, rotating plane of satellites resembling the one in 
Andromeda in cosmological collisionless simulations based on the Cold Dark Matter model. Our new 
high resolution simulations show a correlation between the formation time of the dark matter halo 
and the thickness of the plane of satellites. Our simulations have a high incidence of satellite planes 
as thin, extended, and as rich as the one in Andromeda and with a very coherent kinematic structure 
when we select high concentration/early forming halos. By tracking the formation of the satellites in 
the plane we show that they have been mainly accreted onto the main object along thin dark matter 
filaments at high redshift. Our results show that the presence of a thin, extended, rotating plane 
of satellites is not a challenge for the Cold Dark Matter paradigm, but actually supports one of the 
predictions of this paradigm related to the presence of filaments of dark matter around galaxies at 
high redshift. 

Subject headings: galaxies: individual Andromeda — galaxies: dwarf — galaxies: kinematics and 
dynamics — galaxies: formation — dark matter — methods: numerical 


1. INTRODUCTION 

The success of the currently favored model of struc¬ 
ture formation, the Cold Dark Matter model, lies in its 
outstanding accordance with observations on large scales 
(Tegmark et al. 2004 Springel et al. 2005). However, 
over the last years several problems on galactic and sub 
gala ctic scales have been reported, e.g. the missing satel- 
l ite dKlypin e t al.|1999| iMoore et al.|1999|) and cusp-core 
(Moore 1994 " Boylan-Kolchin et ah] 2011) problems. A 


potential solution to these issues lies m an inappropriate 
comparison between observations and simulations which 
consider only the dark matter component. Cosmological 
simulations including baryonic physics (gas and stars) 
have been shown to be able to ease those problems and 


([Benson et al. [[20071 

Maccib et al.[|2010l Governato et al. 

2010| Di Cintio et al 

..[[2014| Brooks & Zolotov|2014p. 


Andromeda ([McConnachie fc lrwin|2006 Koch & Grebel 


2006 


Ibata eT al.| 2013|), an d possibly 
Way ( jlVletz et al.|2008 ), have sugg 


also around the 

suggested these satel¬ 
lites t o align in a rotationally supported disc. [Ibata 


et al. 


(2013) find that 15 out of 27 satellites in the 


Fan-Andronieda Archa eological Survey (PAndAS, Mc¬ 
Connachie et al. [2009) lie in a thin plane of thickness 
(12.6 ± 0.6) kpc. Furthermore, using line-of-sight veloci¬ 
ties, they report 13 of the satellites in the plane to coro¬ 
tate. This kind of spatial and kinematic align ment is not 
easily found in Cold Dark Matter simulation s ( Pawlowski 
et al. 2014 Pawlowski fc McGaugh|[2014| jlbata et al. 
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2014a|), though no t impossible (Libeskind et al. 2009 


Ibata et al. 2014a). Furthermore there are other stud¬ 


ies which have searched for planes of satell i tes in hydro- 
dyna mical simulations ([Lovell et al.||2011| [Gillet et al. 
2015 and references therein). T'here have been indica- 


tions that filamentary accretion of subh alos can lead to 
anisotropic spatial distributions of them ([Libeskind et al 


2009 Lovell et al. 2011 Libeskind et al. 2014p . JNever- 


theless none of the previous studies was able to closely 
reproduce the parameters of Andromeda’s plane of satel¬ 
lites (number of members, thickness, corotation, and ex¬ 
tension). Since the accretion and distribution of dwarf 
galaxies is primarily governed by the global gravitational 
potential, which is dominated by the dark matter distri¬ 
bution, this issue has raised (again) the question whether 
the Gold Dark Matter model is correct or, conversely, 
needs to be revised. 

In this letter we present new high resolution “zoom-in” 
dark matter only simulations of Andromeda sized dark 
matter halos which reveal thin rotating planes of sub 
halos like in the case of Andromeda. The key insight 
that enables us to solve this puzzle is that halos formed 
through filamentary accretion should form earlier and 
thus can be selected by their higher than average present 
day dark halo concentrations. 

This paper is organized as follows: In §2 we present 
the simulations, including the host halo selection criteria 
and our plane finding algorithm. In section §3 we then 
present the results of our study and in §4 we present our 
conclusions. 


2. SIMULATIONS 

We performed 21 high resolution “zoom-in” dark mat¬ 
ter only simulations using the N-body code PKDGRAV2 
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Fig. 1.— Left panel: Concentration mass rela tion. This diagram show s the concentration as a function of mass of the high resolution 
halos. The solid line is the average relation from|Dutton & Macci6| (|2014|). The dashed line indicates the Icr scatter of this relation. Color 
coding shows the division of the halos into high, average and low concentration. Right panel: Mass growth vs. concentration. The plot 
shows the concentration at z = 0 as a function of mass at 2 : = 2 in terms of the present day mass. Strongly growing halos (late forming) 
are located on the left while least growing halos (early forming) are located on the right. The color coding shows the division of the halos 
into high, average and low concentration. The halos with satellite planes coming closest to the values of Andromeda are marked by black 
circles. The solid one is halo A, the dashed one is halo B. 

equal number of high (red points), average (black points) 
and low (blue points) concentration halos (se e Fig. [^. 

The solid li ne is the power law fit from [Duttcm 
Macci 6 | (|2014), while the dashed ones show the la in- 
trmsic scatter of 0.11 dex around the mean. Our high 
concentration halos have on average an offset of about 2 cr 
from the mean relation. This means these halos are the 
rarest 2.3% of the whole population. In a random sample 
of halos it would thus require ^ 40 simulations to recover 
such rare halo. This helps to explain why previous high 
resolution simulations were unable to reproduce the ob¬ 
served properties of the satellite distribution around the 
Andromeda galaxy: they simply did not sample enough 
halos to find the rarer earliest forming ones. 


( Stadel |2001 2013[ ). We select Andromeda-like mass ha- 
los (5 X 10^^ < M 2 oo/[h~^M q] < 1.5 X 10^^), where the 
halo mass is defined with respect to 200 times the criti¬ 
cal density of the universe. The halos were selected from 
three cosmological boxes of side s 30, 60 and 80 h“^Mpc 
from Dutton & Maccio (2014) which used cos mologi - 
cal parameters from the Planck Collaboration (2014): 

= 0.3175, h = 0.671, ag = 0.8344, n = 0.9624. 

Initial conditions for zoom-in simulations were cre- 
ated using a modified version of the grafic 2 package 


amea version or t ne grahcz packag e 
1) as described in (Penzo et al.||2014|). 
wel was chosen to maintain a roughly 


The refinement level was chosen to maintain a roughly 
constant relative resolution, e.g. ^ 10 ^ dark matter par¬ 
ticles per halo with particle masses of ^ 10 ^ 

This allows us to reliably resolve substructure down to 
^ 5 X lO^h“^M 0 . The force softening of our simulations 
ranges between 0.25 h~^ kpc and 0.36 h~^ kpc. 

2 . 1 . Halo concentration and formation time 

Aside from halo mass, the only other selection criteria 
for our objects is the concentration, which i s a proxy for 
halo formation_time (Wechsler et al. 2002). The right 
panel of Fig. shows the concentration at z = 0 as a 
function of the mass growth since z = 2. The clear corre¬ 
lation validates our approach of using the concentration 
as a first proxy for the halo formation time. 

The reasoning behind such a choice is that, at a fixed 
mass at the present time, early forming halos are more 
likely to form at the nodes of intersection of a few fila¬ 
ments of the cosmic web, while typical halos tend to re¬ 
side inside such filaments (Dekel et al.||2009). One then 
might expect that, rare, early forming halos would ac¬ 
crete satellites from a few streams that are narrow com¬ 
pared to the halo size, while typical halos accrete satel¬ 
lites from a wide angle in a practically spherical pattern. 

In the left panel of Fig. [^we show the concentration- 
mass relation of our high resolution halos. Here the con¬ 
centration is defined as C 200 = 7 ? 2 oo/f- 2 , where R 200 is 
the virial radius, and r _2 is the radius where the logarith¬ 
mic slope of the density profile is — 2 . We select roughly 


2 . 2 . Satellite selection 

Our simulations reveal hundreds of resolved subhalos 
which have to be matched to actual luminous galaxies. 
Galaxy formation models robustly predict the luminous 
subhalos to be the ones most massive at infall times 
(iKravtsov et al.||2004| jConroy et al.||2006| |Vale & Os- 
triker|| 2 U 06 | ). ThW we select a sample of the 3U most 
massive sub halos at the time of the accretion, where we 
restrict our analysis to sub halos within the virial radius 
of the host halo (^ 200 kpc). 

Although observations around Andromeda use a spe¬ 
cial selection function given by the peculiarities o f the 
Pan-Andromeda Archa eological Survey (PAndAS) (Mc- 
Connachie et al. 2009), we choose not to reproduce the 
selection function tor a number of reasons. Firstly, it 
requires surface-brightness information which we do not 
have in our (dark matter only) simulations. Secondly, the 
PAndAS footprint is unique to the Andromeda galaxy, 
being non-circular, and including a region around its 
most massive satellite M33. Thus it would not make 
sense to apply the same footprint to a cosmological sim¬ 
ulation. Thirdly, the spatial depth of the survey is some¬ 
what uncertain due to the difficulty of measuring ac¬ 
curate distances to the satellites. Rather, we apply a 
simple, reproduceable, and physically motivated selec- 



























































Reproducing Andromeda’s Thin Plane of Satellites 


3 


tion criteria. As our satellite population we select the 
most massive sub halos (at the time of infall) within the 
z = 0 virial radius, Rvir- Choosing satellites within the 
virial radius leaves us with a bigger volume Ryir) com¬ 
pared to the observations and hence we use a number of 
30 satellites instead of 27. Furthermore, there is some 
arbitrariness in the number of satellites related to An¬ 
dromeda. Nine known satellites (two that lie inside the 
PAndAS f ield, and seven that l ie outside) were no t con¬ 
sidered by |Ibata et al.| (|2014b| ; jConn et al.| (|2013|) . Nev¬ 
ertheless we experimented witn sample sizes of 25, 27 
and 30 satellites and found no major differences in the 
plane statistics. 


2.3. Plane finding algorithm 

In order to find planes in the distribution of sub ha¬ 
los we generate a random sample of planes defined by 
their normal vector. All planes include the center of the 
main halo. To uniformly cover the whole volume we gen¬ 
erate 100,000 random planes with a fixed thickness of 
2A = 40h~^ kpc. After specifying a plane we calculate 
the distance of every satellite to this plane. A satellite is 
considered to lie in the plane if its distance to the plane is 
smaller than A. For each plane we calculate the number 
of satellites in the plane and its root-mean-square thick¬ 
ness Arms- We then select for every number of satellites 
in the plane the one which is thinnest and richest to an¬ 
alyze for kinematics (for further details see als o |Gillet 
[Ibata et al 


[Conn et ar||2013 ). 


I'he plane of satellites around Andromeda can be char¬ 
acterized by 4 parameters including the number of satel¬ 
lites in the plane (A^in), the number of corotating satel¬ 
lites (A'corot): the thickness of the plane (Arms) and its 
extension. For Andromeda these values are Ain = 15, 
A'corot = 13, Arms = 12.6 ± 0.6 kpc, aud a projected 
diameter of ^ 280 kpc. 


3. RESULTS 

3.1. Plane thiekness vs halo eoneentration 

In a first step of our analysis we investigated the cor¬ 
relation between concentration (as a proxy of halo for¬ 
mation time) and thickness of the plane. Fig. [2 shows 
the plane thickness as a function of the number m satel¬ 
lites in the plane, with lines color coded according to 
halo concentration. There is a clear dependence of plane 
thickness on the concentration of the halo. The thinnest 
planes are only found in the highest concentration (red 
lines), and hence earliest forming halos. Furthermore, 
only high concentration halos have planes as thin as 
observed in Andromeda (assuming 15 members). The 
smooth relation between plane thickness and number 
of satellites in the plane suggests that there is an ar¬ 
bitrariness in the number of satellites chosen to be in the 
plane. We do not see clear evidence of two distinct spa¬ 
tial structures such as a planar and spherical distribution 
of satellites. We further note that an investigation of the 
satellite distribution does not reveal a more concentrated 
satellite distribution in high concentration halos, which 
might trivially explain the dependence of plane thickness 
on concentration that we find. 


3.2. Corotation vs number of satellites in plane 


80 



Fig. 2.— Minimal root-mean-square thickness, min. Arms, of 
planes as a function of number of satellites in the plane. Each 
line represents a different dark matter halo. Red lines show high 
concentration halos, blue lines show low concentration halos and 
blaek lines show average concentration halos. The thinnest planes 
occur in the highest concentration halos. The blue dot shows the 
rms value of the plane of sate l lites o bserved around Andromeda 
(|Ibata et al.|(|2013j);|Conn et al.|(|2013|)). The grey area represents a 
nominal uncertainty m tJie nuniber ot satellites in the plane around 
Andromeda of ±1. 

Fig. 1^ shows the outcome of our plane finding algo¬ 
rithm for the two (early forming) halos best matching 
the values of Andromeda (halo A: Ain = 15, Acorot = 14, 
Arms = 15.8 kpc, diameter ^ 220 kpc; and halo B: 
Ain = 15, Acorot = 13, Arms = 12.9 kpc, ^ 450 kpc). 
The plots show for every value of the number of satel¬ 
lites in the plane, Ain, the number of corotating satellites, 
Acorot- Every dot in the plot represents a different plane. 
In both halos one can find planes with up to 12 members 
that have a 100% corotation fraction. One can always 
find planes with no coherent kinematics (i.e., a corota¬ 
tion fraction of 50%). The points are color coded ac¬ 
cording to the thickness of the plane Arms- As would be 
expected, planes with more satellites tend to be thicker. 
The thickness of the plane is also to first order inde¬ 
pendent of the corotation fraction, except at the highest 
corotation fractions. 

This figure also shows that there is some arbitrariness 
in choosing the best plane from a simulation. For exam^e 
for halo A (left panel in Fig.l^ and upper panel of Fig.®, 
we can find a plane with losatellites and 14 corotatmg 
(one more than Andromeda). If we restrict to a plane of 
14 satellites (instead of 15) we still find a large corotating 
fraction (Acorot = 13) but in a thinner plane Arms = 13.9 
kpc, which is consistent with the Andromeda value. 

We are also able to find a thin plane even richer in 
number of satellites and coherent motion as the one found 
around Andromeda. For example, in halo A we are able 
to find planes consisting of up to 17 satellites from which 
13 share the same rotation direction (see left panel in 
Fig.[^. 

3.3. Visual impression of planes matehing Andromeda 

Fig. i shows one particular projection of two early 
forming halos with plane parameters coming closest to 
the ones observed around Andromeda. The left hand 
side shows the spatial extension and the kinematics of 
the system, clearly revealing a thin but radially extended 
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Fig. 3. — Number of corotating satellites vs. number of satellites in the plane for the two best matching high concentration halos {halo A 
(left) & halo B (right)). The points are color coded by the rms thickness (Arms) of each plane. The square marks the values observed for 
Andromeda (15 in plane, 13 corotating), where the grey shaded area marks the uncertainty of ± 1 satellite in the plane and ± 1 corotating 
satellite. The dots show the values of the planes found for this halo. 


plane with coherent motion of the satellites around their 
host. The right hand side shows a dark matter density 
map of the host halos, with superimposed the two satel¬ 
lite populations, one in the plane (green circles) and one 
outside the plane (black circles). Such di stinct planes 
were reported before in CDM simulations (|Gillet et al 


2015 


Ibata et al.|2Q14a{ Pawlowski et al.|2Q14[ ) but never 


as rich as the ones found here. 

It is worth noticing that previous simulation stud¬ 
ies either were limited by the capability of resolving 
enou gh (sub)structures, e.g., th e Millennium II simula¬ 
tion (Boylan-Kolchin et al.|2009 ) or, conversely, we re lim- 
ited by ver y low number statistics of host halos (Gillet 


et al. 


2015), moreover the halos were not selected accord¬ 


ing to formation time. Since we are able to find planes 
as rich as the one around Andromeda in at least 3 out 
of our 20 simulations, the rarity of the planes can be 
explained by the rareness of early forming halos. 


3.4. Formation scenario 

We now move to the question where does this spa¬ 
tial and kinematic coherence come from. By tracing the 
satellites in the plane back to a redshift of 2 ; = 3 we re¬ 
veal the accretion of the satellites are along dark matter 
filaments. This can be seen in Fig. where we show a 
density plot of the main halo and its substructure at red- 
shift 2 : = 3 indicating the satellites in the plane by green 
circles and the satellites outside the plane by white cir¬ 
cles. Providing two different projections these plots prove 
the connection between accretion along filaments and the 
property of being in a kinematically coherent plane at 
redshift z = 0. 

The presence of a plane of satellites in Andromeda 
seems then to suggest an (unusual) early formation epoch 
for this galaxy. Such a scenario is consistent with other 
observational evidences. For its stellar mass ( lO^^Mo) 


Andromeda lives in a lower mass dark matter halo than 
typical galaxies of the same mass ( Moster et al.|2QI() ). At 
these stellar masses the majority of galaxies are bulge- 
dominated and non star forming, while Andromeda is 
disk-dominated and star forming, consistent with an 
early mass accretion history, devoid of recent major 
mergers. A similar line of reasoning also suggests an 
early formation epoch for the Milky Way. 


4. CONCLUSION 

We have explored the connection between the forma¬ 
tion time of a host dark matter halo and the alignment 
and coherent kinematics of its subhalos using 21 high- 
resolution (10 million particles) cosmological N-body 
simulations. Our key new result is that high concen¬ 
tration (earlier forming) halos tend to have thinner and 
richer planes. Our simulations show that the presence 
of a thin, rotating, and extended plane of satellites like 
the one observed around the Andromeda galaxy is not 
a challenge for the Cold Dark Matter paradigm. Con¬ 
versely it supports one of the key predictions of such a 
model, namely the presence of large filaments of dark 
matter around galaxies at high redshift and the web-like 
nature of cosmic structures in the Universe. 

The connection we report in this letter between the for¬ 
mation time with the satellite distribution at the present 
time, opens a new possibility of constraining the topol¬ 
ogy of the dark matter accretion pattern. 
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Fig. 4.— Edge on view of planes. The left panels show the sign of the line-of-sight velocity of the best matching plane consisting of 15 
satellites (colored dots). Black dots show the satellites not in the plane. The right panels show high resolution density plots of the halo 
hosting the best matching plane. Green eircles indicate the sub halos in the plane and blaek circles indicate sub halos not in the plane. 
The upper panels show halo A, while the lower panels show halo B a slightly more massive halo, also revealing the dependence of radial 
extension on mass. 
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Fig. 5.— High redshift (z = 3) density plots of satellite distribution ending up in the plane and outside of the plane of halo A. The left 
panel shows density plots oi x — y and x — ^-projections of the satellites ending up in the plane at 2 ; = 0 while the right panel shows density 
plots of the X — y and x — ^-projections of the satellites ending up outside of the plane. The upper panel shows that satellites ending up 
in the plane are accreted along two filaments coming from opposite sides of the main halo, which set a preferred infall direction. While 
satellites not ending up in the plane are accreted from everywhere. Comparison with the lower panel shows clearly that satellites ending 
up in the plane lie within the filaments such that their projection collimates in the center of the halo while halos not ending up in the plane 
scatter around the main halo indicating that they are not part of the filaments. 
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